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Abstract 

Dynamics with choice is a generalization of discrete-time dynamics where instead 
of the same evolution operator at every time step there is a choice of operators to 
transform the current state of the system. Many real life processes studied in chem- 
ical physics, engineering, biology and medicine, from autocatalytic reaction systems 
to switched systems to cellular biochemical processes to malaria transmission in 
urban environments, exhibit the properties described by dynamics with choice. We 
study the long-term behavior in dynamics with choice. We prove very general re- 
sults on the existence and properties of global compact attractors in dynamics with 
choice. In addition, we study the dynamics with restricted choice when the allowed 
sequences of operators correspond to subshifts of the full shift. One of practical 
consequences of our results is that when the parameters of a discrete-time system 
are not known exactly and/or are subject to change due to internal instability, or a 
strategy, or Nature's intervention, the long term behavior of the system may not be 
correctly described by a system with "averaged" values for the parameters. There 
may be a Gestalt effect. 

1 Introduction 

Mathematical setting for discrete dynamics is a space X and a map S : X — > X. The 
space X is the state space, the space of all possible states of the system. The map S, 
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the evolution operator, defines the change of states over one time step: x G X at time 
t = evolves into S(x) at t = 1, S(S(x)) at £ = 2, . . . , S n (x) at t = n, etc. If instead of 
one operator, S, we have a choice of evolution operators, S , Si, . . . , Sjv-i> an d at every 
time step choose one of them, then we have a dynamics with choice. One way to visualize 
the multitude of choice through time is to generate the infinite tree of choices. This is an 
infinite rooted tree in which the root has N children, every child has N children, and so 
on. The root corresponds to t — 0, its children correspond to t — 1, the children of the 
children correspond to t — 2, etc. At every step, the children of each node are labeled 
through N — 1. Beginning at the root infinite branches (paths, strategies) represent the 
possible choices: for example, in Figure [T] we choose the path w that starts with Oil... 
(bold edges). For this choice, the first few points in the trajectory of a point x$ G X 
are xi = S (x ), x 2 = Si(xi) = Si(S (x )), x 3 = S 1 (x 2 ) = Si(S 1 (S (x ))), etc. It is 
natural to encode the infinite paths (beginning at the root) by one-sided infinite words 
(strings, sequences) on iV symbols. If w is such sequence, it is convenient to align it with 
the set of non-negative integers Z> and denote by w{k) the {k + l)-st letter of w, i.e., 
w = w(0)w(l)w(2) . . .. Thus, w(0) = 0, w(l) = 1, w{2) = 1, are the first three symbols 
of the path w = Oil . . .. 




Figure 1: The tree of choices in the case of two operators 



• In this paper we study dynamics with choice, i.e., the dynamics of points and subsets of 
X along all possible paths simultaneously. We will explain what this means momentarily. 
Here we would like to emphasize that, from the point of view of long-term behavior, 
dynamics with choice, in general, is not the same as the union of trajectories along different 
infinite paths. We will return to this point later when we talk about the Gestalt effect. 

Let E be the one-sided shift on N symbols, [II]. This means that, first, £ is the set 
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of all one-sided infinite strings w = w(0)w(l)w(2) . . . , where each w(j) is a symbol from 
the list of N symbols [0, 1, . . . , N — 1], and second, there is the shift operator a : E — > E 
acting by erasing the first symbol, a{w) = w(l)w(2)w(3) .... Given the state space X and 
operators So, ... , Sjv-i, we define the corresponding dynamics with choice as the discrete 
dynamics on the state space the product X = X x E with the following evolution operator: 

6: (a^w) i-> (S^o^ar), cr(w)) . (1) 

One can think of a w G E as a plan, a strategy, or as Nature's intervention. Dynamics with 
choice is a language to describe processes where different strategies could be applied or 
happen. Most of mathematical models in natural sciences and engineering are expressed 
in terms of differential equations. Those equations are often continuous limits of discrete 
equations. Continuous case is easier for qualitative analysis. However, there are situations 
where discrete equations describe the processes better. Every realistic model comes with 
parameters. We are interested in situations where parameters may change due to, e.g., 
internal instability or outside intervention. In an illustrative example in section [3j the 
coefficients a and b are proportional to the biting rate of mosquitoes which depends, for 
example, on temperature and humidity which may change from day to day and during 
the day. 

In this paper we study long-term regimes in dynamics with choice. More specifically, 
we define and study global compact attractors in dynamics with choice. By a global 
compact attractor we mean the minimal compact set that attracts all bounded 



sets, see section 2.1 for definitions and references. Thinking in terms of a model with 
parameters, assume we know that for each admissible fixed (in time) set of parameters 
the system possesses a global compact attractor. What happens when the parameters 
switch between admissible values? Is there an attractor? How is it related to attractors 
corresponding to fixed parameters? Is there a Gestalt effect? These are the questions we 
address in this paper. 

There are many real life and engineered systems that switch between different modes 
of operation (the so-called hybrid systems). When the behavior in each mode is modeled 
using continuous dynamics and the transitions are viewed as discrete-time events, such 
systems are called switching or switched. Analysis and especially control of switching 
systems is an area of intensive research, see, e.g., Liberzon's book [20] and the survey by 
Margaliot [22] • There is a natural affinity between switching systems and dynamics with 
choice (see, e.g., [U]), but we will not explore it at this time. 



Readers familiar with iterated function systems, [121 [7J, may wonder if there is a 
connection between iterated function systems and dynamics with choice. Indeed there is, 



but we have to establish it (in section 2.3.3). 
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A (general) Iterated Function System (IFS^can be viewed as a discrete dynamics on 
the space 2 X (of subsets of X). The operators S , Si, . . . , SV-i define evolution on 2 X 
by means of the Hutchinson-Barnsley operator: 

F : A hH. F(A) := S {A) U Si(A) U • • • U S N -i(A) . (2) 

Following a long-standing tradition, people studying dynamics are first of all interested in 
fixed points. In the case of an IFS, those are the fixed points of the Hutchinson-Barnsley 
operator. As has been well illustrated by Barnsley, for many simple IFSs on the plane 
one can (use a computer to) plot their compact fixed points (sets) and obtain fascinating 
fractals, see IE]- Generating fractals is one of the main motivations in the study of 
IFSs. In some papers a fractal is defined as the compact invariant set of an IFS, see [5] 
and references therein. 

To prove that an IFS does have a fixed point, the general definition should be made 
more specific. One needs to specify the properties of the space X; the space 2 X should be 
narrowed to an appropriate class of subsets; assumptions should be made on the operators 
S , Si, . . . , Sjv-i- As an example we state the original result of Hutchinson, [TSJ Section 
3]. 

Theorem (Hutchinson). Let X be a complete metric space (with metric d). Denote by 
B(X) the space of all non-empty closed bounded subsets of X. Assume that each operator 
S , Si, 5tv_i is a strict contraction (i.e., there is a number^ G (0,1) such that 
d(Sj(x), Sj(y)) < jd(x, y) for every pair x,y G X and for all j). Define the evolution 
operator F : B(X) — > B(X) by the formula 

F : A^ F(A) = S {A) U Si{A) U • • • U S N „i(A) . (3) 

Then there exists a unique fixed point K G B{X) of F. Viewed as a subset of X , the set 
K is compact. Also, K attracts every closed bounded subset of X in the sense that, for 
any C G B(X), 

d H {F n {C), K) -> as n -> oo , 
where dn is the Hausdorff distance. 

The IFS with contractive operators are called hyperbolic. Over the years this result 
has been generalized in many different directions (different assumptions on X and/or S^), 
see [1] for references. 

The iterated function systems with probabilities and the 'chaos games' in particular 
show an apparent link to dynamics with choice. Recall that an iterated function system 
with probabilities is an IFS (A; 5*0, ... , Sn-i) together with probabilities p , pi, . . . , Pn-i 

1 We use the abbreviation IFS for single and IFSs for plural forms. 
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assigned to the operators So, Si, ... , Sn-i, where eachp^ > andpo+PH \~Pn-i = 1, 

[7J. The random iteration algorithm (aka the chaos game, [8]) starts with the choice of 
initial state x G X. Next, define recursively x n+ i by choosing its value from the set 
{S (x n ), Si(x n ), . . . , S^-i(x n )} with respective probabilities p , pi, . . . , Pn-i- The choice 
of operators thus will be encoded in some strategy w = w(0)w(l)w(2) ■ ■ ■ G S, i.e., 
x n+ \ = S w f n )(x n ). To show that the sequence (x n ) is determined by xq and w G S, we 
write (x„(x , w)). Consider the averages of the delta-measures concentrated at the points 

x n . It turns out that (under certain conditions) the averages n -1 (5 XQ + S X1 + h 5 Xn _ 1 ) 

converge weakly to the invariant measure of the IFS. More precisely, consider the following 
Markov operator on the space of probability measures on X: 

N-1 

M: v^M(v) = Y,PjS 3 (v), 

3=0 

where Sj(v)(A) = v{Sj 1 {A)) for measurable sets A C X. Hutchinson showed in [13] 
that under the assumptions of the Hutchinson Theorem there exists a unique fixed 
point /i of the Markov operator and the support of /i is the fixed point K of the IFS 
(X; Sq, ... , Sjv-i)- The following theorem was proved by Elton, [10J, with later simplifi- 
cations by Forte and Mendivil, 

Theorem (Elton). Assume X is a compact metric space, the operators Sq, . . . , Sn-i 
are strict contractions on X, and po, Pi, ■ ■ ■ , Pn-i ore the probabilities. Let \x be the 
corresponding invariant probability measure. Then, for any continuous function f : X — > 
E and any xq G X , 



1 n_1 f 

lim - f(xi(x , w)) = / f{x)dn{x) 



for almost all strategies w G S with respect to the product probability measure on S = 
[0, 1, . . . , N — 1] N induced by the distribution [po,Pi, • • • ,Pn-i] on eac h factor. 

Returning to dynamics with choice, we repeat that our interest has not been motivated 
by fractals. We would like to understand the long-term behavior in dynamics with choice. 
We assume that X is a complete metric space (with metric d), the operators Sq, . . . , S^r-i 
are continuous, and each of the (semi) dynamical systems (X, d, Sj) possesses a global 
compact attractor. Consider the corresponding dynamics with choice as the dynamics on 
the product metric space |l = IxE generated by the operator (5 acting according 
to the rule 0. From general theory (see section 2.1) we know that a system ought to 



2 Yi can be equipped with a metric making it a compact metric space, see Section 2.3.1 for a specific 
choice. We denote here by dist the corresponding product-metric on X x E. 
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enjoy certain compactness and dissipativity properties in order for it to possess the global 
compact attractor. 

In general, even when the individual systems (X, d, Sj) do have attractors, the system 
(X, dist, &) will not have a global compact attractor. There are several reasons why. One 
counter-example we borrow from [3] (where it is used in the context of IFS). Take X — R 
with standard metric d and define two maps, So and Si, as follows: 



Each of the systems (X, d, Sj) has the global compact attractor, a singleton {0}. At the 
same time, the trajectory x n = S^n-i) o 

Sw{n-2) o • • • o S w (o)(xq) corresponding to the 
periodic string w = 010101 ... is unbounded for any initial point xq ^ 0. Hence, there is 
no compact attractor attracting (xq,w). 

The second example is infinite-dimensional. Let Bq = Bq(pq) an d £>i = B\(pi) be two 
disjoint closed unit balls centered at po and pi in an infinite-dimensional Banach space. 
Let X = B Q UBi. Define the maps So and Si as follows: on B the map Sq is a contraction 
and it maps B\ to B ; the map Si is a contraction on B\ and maps B to B\. 



The system (X, d, So) does have the global compact attractor, {po}, and (X,d,S\) does 
have the global compact attractor, {p±}. The corresponding dynamics with choice, 
(X, dist, &), does have the global closed attractor, namely, X, but does not have the 
global compact attractor. 

In the first example, the maps are compact (which is good), but they do not have a 
joint bounded absorbing set (lack of dissipativity in (X, dist, &)). In the second example, 
there is a joint bounded absorbing set, B U B 1: but there is not enough compactness (the 
maps Sj are not compact, not contracting, and, more generally, not condensing). 

These examples show what kind of situations do not allow global compact attractors 
in the dynamics with choice. Thus, we make additional assumptions. First, we assume 
that there exists a bounded absorbing set that absorbs every bounded set regardless of 
the strategy. In applications, absorbing set is usually a ball of the radius that depends 
on the parameters of the model. Our "dissipativity" assumption means that there is a 
common estimate on the radius for different values of the parameters. 

Assumption 1. There is a closed, bounded set B C X such that for every bounded A C X 
there exists m(A) > such that S w / n -i) ° 5 , u ,(„„ 2 ) o • • • o S w ^ (A) C B for every word 
w = w(0)w(l) . . . w(n — 1) of length n > m(A). 
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Our second, "compactness" assumption is that each of the operators Sj is condensing 
with respect to a common measure of noncompactness. This assumption covers practically 
all situations encountered in applications: contractions, compact operators, and compact 
plus contractions. As their name suggests, measures of noncompactness measure how far 
a set is from being compact. There are several different measures of noncompactness in 
use, pp. For example, the Hausdorff measure of noncompactness of a set A is the infimum 
of e > such that A has a finite e-net. In this paper we use only very general properties 
shared by all popular measures of noncompactness, see Definition [7] in section 2J2_ below. 



Let ip be a measure of noncompactness (as in Definition [7]). An operator 5* : X — > X 
is condensing with respect to if) iff if>(S(A)) < if) (A) for any non-compact set A, and 
if)(S(A)) = if) (A) = if A is compact. Our second general assumption is this. 

Assumption 2. Each operator Sj is if) -condensing. 



In section 2.3 we prove the following result. 



Theorem 1. Let X be a complete metric space and let So, Si, . . . , Sn-i be continuous, 
bounded (i.e., take bounded sets to bounded sets) maps X — > X . In addition, let assump- 
tions 1 and 2 be satisfied. Then the system (X, dist, &) has a global compact attractor, 
97t. (That 971 is the global compact attractor means that 971 is the smallest compact in X 
attracting every bounded set in X.) 

The attractor 971 has the following properties. 

(1) 971 is (strictly) invariant: (5(971) = 97T. 

(2) 97t is the union of all closed bounded sets A G X with the property A G &{A). 

(3) 971 is the maximal closed set with the property A C &(A); in particular, 971 is the 
maximal (strictly) invariant closed set. 

(4) Through every point (x,w) G 97T passes a complete trajectory. This means there ex- 
ists a two-sided sequence . . . , X-2, X-i, Xq, xi, X2, ■ ■ ■ of points in X and a two-sided 
infinite string . . . s{— 2)s(— l)s(0)s(l)s(2) . . . such that x(0) = x and s(0)s(l)s(2) • ■ ■ = 
w(0)w(l)w(2) . . . and such that S s ^(x n ) = x n+ i for every integer n. 

(5) 971 is the union of all complete, bounded trajectories in X. 

Given the state space X and operators So, ... , SV-i, there are two ways of describing 
dynamics generated by the corresponding IFS. First, one can follow the trajectories of 
bounded subsets of X under the iterations of the Hutchinson-Barnsley map F, see 
We denote such system by (X, d, F) . The notion of the global compact attractor as the 
minimal compact set that attracts all bounded sets, is well-defined for (X,d,F). The 
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second possibility is to choose the space of closed bounded sets, B(X), as the state space 
of the system and study the dynamics of its points under the iterations of F. As a rule, 
B(X) is equipped with the Hausdorff distance cf#. Thus we obtain the second system, 
(&(X),cLh, F). It turns out that from the point of view of global compact attractors 
the dynamical system (B(X),dn, F) is not very interesting (because convergence in the 
Hausdorff metric is too strong). It possesses an attractor (in the sense we use here) 
essentially only if the maps Sj are contractions, so then the attractor is just one point in 
B(X). For more general Sj, it makes more sense to study the fixed points of F. 

In sections 2.3.4 and |2.3.5 we establish the following connection between the dynamics 



with choice and the corresponding IFS. 

Theorem 2. Make the same assumptions on the space X and operators So, ... , SV-i as 
in Theorem^ Then 

(1) The IFS (X,d,F) does have a global compact attractor, K. 

(2) The set K is the largest compact set in X which is invariant under the Hutchinson- 
Barnsley map F , K — F(K). 

(3) The attractor DJt of the dynamics with choice has the following product structure: 

m = k x s. 

In the extensive literature on IFSs the main question is the existence of "the fractal", 
i.e. the maximal compact set invariant under the Hutchinson-Barnsley operator F. This 
corresponds to the second assertion of our Theorem [2] We believe that viewing "the 
fractal" of an IFS as the attractor of the dynamical system (X, d, F) is beneficial to the 
theory of IFSs. This approach, in particular, points to the "right" assumptions on the 
space X and the operators Sj. 

Iterated function systems with compact (possibly multi-valued) operators have been 
considered previously, see, e.g., [3]. The statement of Theorem 5.8 in [3J which establishes 
the existence of a compact set invariant under F, needs some additional (dissipativity) 
assumption such as our Assumption 1, for example. The IFS with condensing (and multi- 
valued, in addition) operators have been considered by Lesniak, [19] and Andres et al., 
[3]. The assumptions of Theorem 3 in [4] require that the image F(X) of the whole space 
X be bounded. The word "minimal" referring to the "fractal" in [U Theorem 3] should 
probably be replaced by "maximal", see also |19j Theorem 3]. 

Our assumptions on the state space and the operators guarantee that, for every fixed 
j = 0, . . . , N — 1, the discrete dynamics generated on X by Sj does possess the global 
compact attractor (in X). More generally, as we show in sections 2.3.4 and 2.3.5, it 



makes sense to define individual attractors, A w corresponding to every string (infinite 
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path in the tree of choices) w G E. The attractors generated by each Sj correspond 
to "constant" strings, w = jjj .... It is not hard to see that such attractors do not 
exhaust the attractor (fractal) K. There are situations when the union of all A w is K 
(this happens, in particular, when S^s are strict contractions). However, in general, the 
union |J A w is strictly smaller than K. We give an example of this in section 12.3.5 In 



the cases when |J A w is strictly smaller than K we say that there is a Gestalt effect, 

i.e., "the whole is greater than the sum of its parts." This is a new phenomenon. It has 
not been observed in the framework of Iterated Function Systems because, as we show in 



Lemma 15 , the Gestalt effect cannot occur when operators Sj are contractions. 



An important generalization of dynamics with choice is dynamics with restricted choice. 
The name should indicate that not all strategies (sequences w = w(0)w(l) ■ ■ ■ G E) are 
allowed. In particular, we consider the sets in E that are closed and shift invariant, 

1. e., subshifts, see j2Tj [13]. Given a subshift A C E, we consider the dynamics on the 
product-space Xa = X x A generated by the map & as in Q. 

Theorem 3. Let the space X and the operators So,- ■ ■ , SV-i satisfy Assumptions 1 and 

2. Let A be a (one-sided) subshift o/E. Consider the dynamical system (Xa, dist, &). 

(1) The dynamical system (Xa, dist, &) does possess a global compact attractor, 971a. 

(2) The attractor 971a is invariant in the sense that (5(97Ta) = 971a- In fact, 971a is the 
maximal invariant compact set in Xa- Also, 971a is an invariant compact subset of 
the global attractor 971 of the unrestricted dynamics (X, dist, &). 

(3) Through every point (x(0),w) G 97Ta passes a complete trajectory, i.e., there exist 
a two-sided sequence of points . . . , x(— 1), x(0), x(l), . . . and a two-sided symbolic 
sequence . . . w(—l)w(0)w(l) . . . extending w (in the extension of the subshift A) 
such that S w (i)(x(i)) = x{i + 1) for all integers i. 

(4) Let K\ denote the projection of the attractor 97Ia onto the X component. The 
set K\ is a compact subset of the set K of Theorem [S| There exist compact sets 
A , . . . , A N _i such that K\ = A U A\ U ■ ■ ■ U A N _i and 

K A = A U A 1 U ■ • ■ U A N ^ = S (A ) U S 1 {A 1 ) U • ■ ■ U S'at-i^-i) (4) 

(5) In general, the attractor 971a is not a product. There may be infinitely many different 
sets among the slices M.{w) = {i6l : (x, w) G 971a}- However, if the subshift A 
is sofic, the number of different slices is finite. 
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This theorem is proved in section 2.4 For information on sofic shifts see [2T] and section 

Restricted dynamics of a sort has been considered previously, see [2UI23]- For example, 
the graph directed Markov systems of [23] describe iterations of uniformly contracting 
maps indexed by the edges of a directed (possibly infinite) graph. In this case there is a 
correspondence between the points of the limit set and the infinite walks through the graph 
(the coding space). Similarly, the directed IFSs discussed in [8] are defined with the help 
of the aforementioned correspondence, and the fractal (or attractor) is understood in 
terms of the map from the code space to K as the image of A, [SJ Theorem 4.16.3]. The 
correspondence between the points of the code space, E, and the points of K is possible 
because the maps are contractions (right away, or eventually). Our approach gives a new 
and more general view on restricted dynamics. We justify the name - attractor - and 
unveil attractors' more subtle structure (assertion 5). This new approach allows us to 
work in a much more general setting and with transformations that are not contractions. 
We do not have and do not use a map from the code space into the attractor. 

We should mention the paper of Andres and Fiser, [2]. They use their result of [3] on 
the existence of the fractal (the set K in our notation) for an IFS with compact opera- 
tors Sj to conclude that fixed time solution operators of systems of ordinary differential 
equations could play the role of maps generating the IFSs. As an illustration they use 
five two-dimensional systems of ODEs to produce five operators (incidentally, contrac- 
tions, as noted in [2]) and plot the corresponding dragon-tail-like fractal set. Although 
their message is that IFSs and fractals can be generated by solution operators of ODEs, 
their examples can serve as an illustration for our dynamics with choice attractors (due 
to Theorem 2(3)). 




Figure 2: Attractors for (X, d, S ) (right) and (X, d, Si) (left). 

In section [3] we apply the theory to a specific example of a discrete Ross-Macdonald 
type model of malaria transmission. The model can be viewed as a time discretization 
(with time step At) of the ODE model, or as a pre-ODE form of the model. The reason 
we have chosen this model is because it is simple and we can visualize all the attractors. 
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The state space is the unit square X = {(x,y) : < x < 1, < y < 1}. We use 
two sets of parameters, which define two operators, So and Si. Those operators are not 
contractions, but they are compact, because the system is finite-dimensional. The discrete 
dynamical system generated on X by S has two fixed points, (0,0), which is unstable, 
and (11/15, 11/16), which is stable. The system generated by 5*1 also has two fixed points, 
(0,0) (unstable) and (7/25, 7/12) (stable). The attractors of the systems (X, d, So) and 
(X, d, Si) are just the heteroclinic trajectories connecting the unstable and stable fixed 
points. They are depicted on figure [2j The effects of freedom of choice on the dynamics 
are as follows. For the unrestricted dynamics with choice, we see (figures E] and [BT) that 





Figure 3: 
At = 0.05. 



The attractor slice K: 



Figure 4: The attractor slice K; 
At = 0.005. 



the attractor, K, is a rather big set with the two individual attractors corresponding to 
So and Si, respectively, forming parts of the boundary of K (the right and left sides). The 
remaining part of the boundary is quite irregular when At is relatively large, figure [3j For 
smaller At, this part of the boundary becomes much smoother and looks like a smooth 
curve, figure |4j In the limit At — > 0, the set K retains its two-dimensional fullness. It 
is not an attractor of an ODE with averaged parameters. In fact, all such attractors are 
one-dimensional (each is a heteroclinic trajectory connecting two fixed points) and lie 
inside of K. 

We consider also the dynamics with restricted choice corresponding to the golden mean 
subshift. We exhibit two different slices in the global attractor. Their projections on the 
state space do overlap, and their union is smaller than the set K for the full shift. 

Finally, we remark that more general compact and condensing operators will be needed 
in the study of dynamics with choice related to nonlinear dissipative partial differential 
equations, which we plan to address at a later time. 

The structure of the paper. This Introduction is followed by two chapters. In 
the long chapter [2] we address theoretical questions. In section 2.1 we give definitions 



and state the basic results related to global compact attractors. Section 2.2 deals with 



measures of noncompactness. In section 2.3 we study dynamics with choice. Dynamics 



11 



with restricted choice is studied in section |2.4 In chapter [3] we analyze a simple example 
illustrating some of our theoretical results. Despite its simplicity, this example shows that 
the result of dynamics with choice is larger than the sum of its parts. 



2 Theory 

2.1 Attractors: general facts 

We start by collecting the basic facts about attractors. There are several books such as 
H21 El HS1 EE] devoted to this subject. Our presentation is closer to [T7]. We present 
only the results that we need. For the proofs of Theorems [5] and [6] see the books quoted 
above. 

Let Y be a complete metric space with metric dist and let $ : Y — > Y be a continuous 
map. Iterations, $ n , of $ define a discrete (semi) dynamical system on (Y, dist). It is useful 
to consider not only the dynamics of individual points under the action of $, but, more 
generally, the dynamics of bounded sets. Denote by B(Y) the collection of all bounded 
subsets of Y. We say that the set A G BiY) attracts the set B e B(Y) if 

dik ($"(£), A) -> 0, 

n— >oo 

where the one-sided distance between two sets, dist (C, A), is understood as sup dist (y, A). 

Definition 4. We call a set 971 C Y the global compact attractor of the system (Y, dist, $) 
if 

• 371 is compact, 

• 971 attracts every bounded subset of Y , 

• 971 is the minimal set with these two properties. 

For a system to possess a global compact attractor, it should enjoy certain properties, 
namely, some form of compactness and some dissipativity. Here is the basic existence 
(and uniqueness) result. 

Theorem 5. The semidynamical system (Y, dist, $) has a global compact attractor if and 
only if it enjoys the following two properties: 

1. ("compactness") For every bounded sequence {y^) in Y and every increasing sequence 
of integers n k — > +oo, the sequence $ nfe (|/fc) has a convergent subsequence. 
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2. ("dissipativity") There exists a bounded set BcF which absorbs every bounded set 
in the sense that for every A e B(Y) there exists m(A) > such that $ n (A) C B 
for all n > m(A). 

Some of basic properties of a global compact attractor are collected in the following 
theorem. 

Theorem 6. Assume that fXft is the global compact attractor of the semidynamical system 
(Y,dist,<$>). Then 

1. 9Jt is the union of all possible limits of sequences of the form $ nk (yk), where yj. is a 
bounded sequence in Y and nk — > oo. 

2. 9JI is (strictly) invariant: $(£DT) = DJl. 

3. 9Jt is the union of all closed bounded sets A with the property A C <&{A). 

4- 3JI is the maximal closed set with the property A C in particular, DJl is the 

maximal (strictly) invariant closed set. 

5. Through every pointy 6 DJl passes a complete trajectory, i.e., there exists a two-sided 
sequence . . . , y_ 2 , y~i, yo, Vi, ■ ■ ■ of points in DJl such that y = y and y m +i = &(y m ) 
for all integers m. 

6. OJl is the union of all complete, bounded trajectories in Y . 

In applications, people do not verify the "compactness" property of Theorem [5] directly. 
Instead, they use one of the known sufficient conditions that imply it. Two of the most 
useful sufficient conditions are: 

• $ is a compact map (i.e., $ : Y — > Y is continuous and maps bounded sets into 
relatively compact sets); 

• in the case Y is a Banach space, $ is a sum of a compact operator and a strict 
contraction. 

Compact $ arise, e.g., in the finite-dimensional dynamics described by differential or dif- 
ference equations, or, in the infinite dimensional case, in dynamics described by parabolic 
equations. The "compact + contraction" $ appear, e.g., in hyperbolic problems with 
damping. Each of the two sufficient conditions implies that $ is condensing with respect 
to some measure (s) of noncompactness. Since measures of non-compactness and condens- 
ing operators are not widely known, below we give a brief account of the facts we need 
and refer to pQ for more details. 
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2.2 Measures of noncompactness 

Measures of noncompactness assign real non-negative numbers to bounded sets with value 
assigned exclusively to relatively compact sets. The basic examples are the Kuratowski 
measure of noncompactness a and the Hausdorff measure of noncompactness \. By 
definition, a (A) is the infimum of numbers e > such that A admits a finite cover by 
sets of diameter less than e. The number x(A) is the infimum of those e > for which A 
possesses a finite e-net in Y. In this paper we adopt the following definition of a general 
measure of noncompactness (our definition differs from that in [TJ). 

Definition 7. A function ip assigning non-negative real numbers to bounded subsets of 
(a complete metric) space Y will be called a measure of noncompactness iff it has the 
following properties: 

(i) ip(A) = if and only if A is relatively compact; 

(ii) If A x c A 2 , then ^(Ai) < if>(A 2 ) ; 
(%u) %l)(Ai U A 2 ) = max {V>(^i), if>(A 2 )} / 

(iv) There exists a constant c(ip) > such that 

|^(A)-^(A 2 )| < c^)d H {A x ,A 2 ), 
where dn is the Hausdorff distance, 

d H (A h A 2 ) = max {dist (A u A 2 ), dist (A 2 , At)} . 



Both a and x en j°y all these properties. Note that property (iv) implies that the measures 
of noncompactness of a bounded set and its closure are equal: 

(v) V(3)=V04). 

Definition 8. A continuous bounded map $ : Y — > Y is called condensing with respect 
to the measure of noncompactness tp (we also say $ is ip -condensing) iff ip($(A)) < ip(A) 
for any bounded A, and ip(&(A)) < ip{A) if ip(A) > (i.e., if A is not compact). 

Theorem 9. Consider the system (Y, dist, $). Assume that $ is condensing with respect 
to some measure of noncompactness ip and that there exists a bounded set B which absorbs 
every bounded set. Then (Y, dist, $) possesses a global compact attractor. 

In the case ip is the Kuratowski measure of noncompactness this result is proved by 



Seda, [25]. In Lemma 10 below we establish a more general result. 
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2.3 Dynamics with choice 

2.3.1 Words, strings 

Fix an integer N > 1. Using the integers 0, . . . , N — 1 as the alphabet, construct strings 
(words) of finite length and (one-sided) strings of infinite length. Denote by E* the set 
of all finite length strings (words), and denote by E the set of all (one-sided) infinite 
strings. The word of length is the empty word. The set of non-empty words is denoted 
E + . Given a string w G E* U E, w(0) is the first letter of w, and w{k) is the {k + l)-st 
letter of w. The length of w is denoted \w\. If w is a finite string and u G E* U E, their 

concatenation is denoted w.u; if \w\ = m, then {w.u){m + k) = u(k) for k — 0, 1, For 

a u> G E* and s G E* U E, we write u> C s if w is the beginning of the string s, i.e., if there 
exists u G E* U E such that s = w.u. For an infinite string, s, its first n letters form a 
word denoted s[n], i.e., s[n] = s(0)s(l) . . . s(n — 1). The set of all words of length m will 
be denoted by E^. 

Equip the space E*UE with the metric where dj:(u, v) = 2~ m if u[m — 1] = v [m — 1] 
but u[m] 7^ v[m\. It is well-known, [8], that both E* U E and E with metric ds are 
compact. The shift operator, a, acts on infinite strings by deleting the first letter, i.e., 
<y(u) = u(l)u(2) .... The shift operator maps E onto itself. It is continuous; in fact, 
d s (a(u), a(v)) < 2d s (u, v). 

2.3.2 The skew-product dynamics 

Let X be a complete metric space with metric d, and let So, Si, ... , Sjv-i be continuous, 
bounded maps X — > X. Define the product metric space X = X x E with metric dist, 

dist ((x, u), (y, v)) = d (x, y) + d s (u, v) . 

The skew-product dynamics on X is generated by the map & : X — ► X acting according 
to the rule 

& (x,u) = (S^o) (x), cr(u)) . 

This map is obviously continuous and bounded. Because we will consider iterations of & 
such as 

& l (x, U) = (S^n-l) O • • • O O 5 U ( ) (x), (T n (u)) , 

we introduce the notation 

S w = S^n-i) ° • • • ° 'S'tu(i) ° ^(o) j 

where w is a word of length n. Thus, we can write & n (x, u) = (S u [ n ] (x), a n (u)). 

Assumption 1. Assume there is a closed, bounded set Bel such that for every bounded 
A C X there exists m(A) > such that S w (A) C B for every word w of length n > m(A) . 
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[In applications B is usually a closed ball of radius that depends on the parameters of 
the model. Showing that for different values of the parameters there is a common estimate 
on the radius is enough to verify Assumption 1.] 

Let if) be a measure of noncompactness as in Definition [7j 

Assumption 2. Assume that each operator Sj is ip- condensing. 

We are going to apply Theorem [5] to prove the existence of the attractor. For this to 
work we need to justify the following fact: 

For every bounded sequence (xk,Uk) G X and every increasing sequence of integers 
n k — > +oo, the sequence & nk (xk,Uk) has a convergent subsequence. 

Thus, pick a bounded sequence G X and any sequence u k G S. Pick an increasing to 
+oo sequence 7i&. Since £ is compact, the sequence a nk {uk) has a convergent subsequence. 
So, we may assume from the very beginning that a nk {uk) — > s G S. Denote Wk = Uk[rik\- 
We have 

e nk (x k ,Uk) = (S Wk (x fc ), o nk (u k )). 

Since the E-component converges, we need to show that the sequence S Wk (x^ in X 
has a convergent subsequence. Because £* U S is compact, we can choose a convergent 
subsequence from Wk- We will assume that Wk itself converges to some w6E. 

Lemma 10. Under the Assumptions 1 and 2, let w n be a sequence of finite words of 
increasing lengths \w n \ — > +oo and w n — > w G S. For any bounded sequence x n the 
sequence S Wn ( x n) has a convergent subsequence. 

Proof. Pick a bounded sequence Xk- By Assumption 1, when the length of the word 
w n is sufficiently large, S Wn ({xk}) C B. Dropping the first few terms if necessary, we 
will assume that S Wn ({xk}) C B for all n. Next, since w n — > w, for arbitrarily large m* 
we have io n [mJ = io[m*] for all sufficiently large n. Choosing m* large enough, we have 
S w [m,]{{xk}) C B. Writing w n = w[mJ.s n , we see that 

Sw n ({Xk}) = S Sn (S w[m j({x k })) 

and s n — > s = cr m *(u;). Use Assumption 1 to find m(B) and assemble the set 

B= |J 5„(B). 

The set B has the property that 5^(B) C B for all finite words v. If we choose m* above 
sufficiently large (first choose it to guarantee S w [ mif \ ({<£&}) C B and then increase it by 
m(B)), the set S^m*] ({xfc}) will be inside of B. We will reformulate our problem now. 
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We have a sequence yu = ^[m.]^*) i n B and a sequence s& — >■ s. We need to show that 
the sequence S Sk (yk) is relatively compact. 

Consider the positive (infinite) trajectory of B under the Hutchinson-Barnsley evolu- 
tion: 

C = BU |J |J S V (B). 

Define inductively 

j=0,...,N-l 

We have C D C\ D C 2 D 

Introduce a collection Sj of all sets icB that can be represented in the form 

A = A n , where A n is a /mrfe (or empty) subset of C n . 

n>0 

We show next that every set A £ S) is relatively compact. Because the sequence {S Un (y n )} 
is in fj, our lemma will be proved. The argument that follows is a modification of a part 
of [U Lemma 1.6.11]. 

First, we claim that there exists a set A* £ such that 

ijj(A*) = sup ^(A) • 

This follows from Lemma 1.6.10 of pp. Although their lemma is stated for the Hausdorff 
measure of noncompactness, their proof uses only the properties (i), (ii), and (iii) of 
Definition [7] and, therefore, works for our if>. 

Now, A* = |J A* n1 where each A* n is a finite (or empty) subset of C n . For every 

n>0 

pG A* with n > 1 pick a g 6 C n -\ and a j £ [0, 1, . . . , N — 1] such that Sj(q) = p. Denote 
the resulting subset of C„_i by A n _i and define A = [j A n . Clearly, A £ Sj. Hence, 

n>0 

ip(A) < if>(A*). Now, consider the set F(A) = S (A) U ■ ■ ■ U S N -i(A). Clearly, F(A) £ £ 
and F(A) D U n >i Hence, 

V>(A*) = V(U A;) < ^(F(A)) = maxMSo(A)), . . . , ^-i^))} < . 

n>l 

The properties of if> used here are, from left to right: the first equality uses properties (i) 
and (iii), the first inequality follows from (ii), the second equality follows from (iii), the 
second inequality is due to the assumption that each operator Sj is ^-condensing. Now, if 
A is not relatively compact, we must have ifj(Sj(A)) < if} (A) for every j. This would imply 
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if) (A*) < if){A) which contradicts the fact that if) (A) < if>(A*). Thus, if){A*) = if)(A) = 0, 
and so every set in .f) is relatively compact. This concludes the proof of Lemma. 

Applying now Theorems [5] and [6j we immediately obtain Theorem [TJ We next proceed 
to the proof of Theorem [2] 



2.3.3 Proof of Theorem H 

Consider the IFS dynamics (X,d,F). Having made Assumption 1 we will prove the 
existence of a global compact attractor K C X if we show that for every bounded sequence 
(x k ) C X and every increasing sequence n k — > +oo, the sequence F nk (x k ) is relatively 



compact. But this follows from Lemma 10 Thus, the existence of the global compact 



attractor for the IFS is established. The global compact attractor, K, comes with all the 
properties listed in Theorem |6j In particular, K is the maximal compact set invariant 
under F. 

To prove that 971 = K x X we start by showing that the slices of the attractor corre- 
sponding to different strings are all the same, i.e., the set {x G X : (x, s) G 971} does not 
depend on s. 

All slices are equal. Recall that every point (x, s) in OK is a limit of some sequence 
& nk (x k , s k ) with bounded (x k ) C X and cr nk (s k ) converging to s. As we argued above, 
we can write 

e nk (x k ,s k ) = (S Wk (x k ),a nk (s k )), 

where w k is a prefix of length \w k \ = n k of the string s k , i.e., s k = w k .a nk (s k ). The 
sequence S Wk (x k ) converges to x and <J nk (s k ) converges to s. The limit of the pair will 
not change if we replace s k by w k .s. Clearly, for any string u G S, we have 

lim (S Wk (x k ),a nk (w k .u)) = (x,u) . 

This proves that DJl = A x S, with a compact set A C X. 

Since S = 0.SU1.EU- • -U(JV-l). Sand since 6(071) = 07t, we get &(Ax E) = (S (A)U 
S^A) U • • ■ U S N -i(A)) x S = A x E. In other words, A = S (A) USi(A) U • • • U S N -i(A). 
Because K is the maximal compact in X with this property, we have A C K. On the 
other hand, &(K x S) = K x E. Since A x E is the maximal compact in X with this 
property, we have K C A, and hence, A = K. This completes the proof of Theorem [2] 

2.3.4 Individual attractors 

Every fixed strategy also generates a dynamics on X: if w G E is the (fixed) strategy, then 
an x G X moves to S w (p)(x), then to S^i) (S w ( )(x)) , then to S w ( 2 ) (-S'w(i) {S w (o){x))), etc. 
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Denote this dynamics by (X,d,w). This is not a (semi) dynamical system, but we should 
not worry about names. Certain important notions related to the long-term behavior with 
natural adjustments still make sense. For example, the individual, i.e., corresponding to 
an individual strategy w, trajectory of a set B is the union 

BUS w[1] (B)US w[2] (B)U... . 

We define the individual u-limit set of a bounded set B as 

u(B, w) = {y EX : y = lim S w[nk] (y k ) for some sequence (y k ) in B} . 

By analogy with Definition |4| we say that a set A is the global compact attractor of 
system (X, d, w) if it is the minimal set with the following two properties: A is compact 
and A attracts every bounded set under the strategy w, i.e., for any bounded B, we have 
lim^oo dist (S w [ n ](B), A) = 0. 

Next theorem establishes the existence of individual compact attractors, A w , of sys- 
tems (X, d, w). Along the way we establish various properties of the u;-limiting sets 
uj(B, w). 

Theorem 11. Under the Assumptions 1 and 2, every system (X,d,w) has the global 
compact attractor, which we denote by A w . This attractor is the intersection of the closures 
of the tails of the trajectory of the absorbing set B ; 



•Aw — Pi [J S w [ k ](B) 



n>l k>n 

The attractor, A w , is the union of all u(B,w) with bounded B. 

Proof. We use some notation and keep in mind the argument from the proof of Lemma 



10 Due to Assumption 1, every bounded set eventually finds itself in the set B and after 
that stays there. 

Step 1. The u;-limit sets of bounded sets are not empty. 

Pick a point x G X and follow its trajectory, x n = S w [ n ](x ). There will be a time n 
such that x n G B C C , and then inevitably x n+1 G Ci, x n+2 G C 2 , and so on. By 



Lemma 10, the sequence (x n ) is relatively compact. Thus, uj({x q },w) ^ 0. Because 
uj({xo},w) C u(B,w) if xo G B, we have u(B,w) ^ 0. 

Step 2. u(B,w) is the intersection of the closures of the tails of its trajectory, 
hence cu(B,w) is closed. 

Note that u(B,w) can be characterized as follows. uj(B, w) is the set of all y G X such 
that for every e > and every integer k > there exist and x G B and n > k so that 
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S w [ n ](x) E O e (y) (where O e (y) is the e-neighborhood of y). Yet another way to describe 
u>(B, w) is to consider the trajectory of B and its tails: 

D = B U S w[1] (B) U S w[2] (B) U . . . , D n = (J (£) . 

m>n 

Clearly, D D Di D D 2 D It turns out that 

u{B, w) = f| W n . (5) 

n>0 

Indeed, inclusion C is obvious. To prove the "D" part, pick a y in the intersection of 
the tails and set e n = 2~ n . In Di there is a point y\ = S w [ mi ](xi), x\ E B, such that 
d(yi, y) < ei. In .D mi +i there is a point yi = S w [ m2 ](x 2 ), rn 2 > mi, such that d (y2, y) < £2, 
and so on. The limit of y n belongs to u(B, w), i.e., y E u>(B, w). 

Step 3. Lo(B,w) is compact. 

Compactness of u(B, w) will follow from the fact that the intersection of the closures 
of the sets C n in the proof of Lemma [10] is compact, because, thanks to Assumption 
1, P| D n C P| C n . Denote C* = f] C n . Since u(B, w) is not empty, C* is not 

n>0 n>0 n>0 

empty as well. And it is closed. If the set C* is not compact, then there exist e > 
and an infinite sequence (y n ) C such that d(y n , y m ) > e for all n and m ^ n. 
Since y n E C n (in fact, the whole sequence lies in every set C n ), there exists a sequence 
y n k E C n that converges to y n as k — >■ 00. For every e > there are numbers k n such that 
^ (?/nfc„? Vn) < e f° r ah n - When e < e /2, the Hausdorff distance between the sets {y n k„} 
and {y n } is n ot greater than e. Using property (iv) of the measure of noncompactness ijj, 



we obtain \ip ({y n }) - $ ({y n k n }) I < c(^) e. Now, ^ ({ynfc„}) = by Lemma [10| Then 
^ ({Vn}) < c(-0) e. Since this is true for any e, we obtain ip ({y n }) = 0, a contradiction. 
This proves that is compact. 

Step 4. uj{B, w) attracts B. 

To show that for every e > there exists an m such that S w [ n ](B) C O e (cu(B, w)) for all 
n > m, we argue by contradiction. Assume there exists an eo > such that S w [ n ](B) does 
not lie inside O tQ (u)(B, w)) for infinitely many n. This means that there is a sequence 
Xk in B and a sequence — > +00 such that 5 , w [ nfc ](xfc) ^ O eo (u(B, w)). But we already 
know that S^n*.]^) must have a convergent subsequence whose limit must be in u(B, w). 
A contradiction. 

Step 5. A w = u>(B, w) = |J u{B,w). 

bounded B 

Because every bounded set is eventually absorbed by the set B, we have u>(B, w) C 
uj(B, w). Thus, uj(B, w) attracts every bounded set. It is compact and minimal, hence, 
it is the global compact attractor. The theorem is proved. 
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2.3.5 Interplay between individual attractors 

Recall, that (with Assumptions 1 and 2) the global attractor DJl of (X, dist, E) is a product 
971 = K x E. 

We start with a few simple observations. 

Lemma 12. A w C F (A w ) C K , where F is the Hutchinson-Barnsley operator. 

Proof. Pick a point, x, in A w . Then x = lim S w [ nk ](xk) for some bounded sequence (xk) 
in X and — ► oo. Among the last letters of the words w[rik] there is at least one, say, j, 
that repeats infinitely many times. Sparse the sequence so that every w[n k ] has the last 
letter j. Then, 

S w [n k ]( x k) — Sj (S w [ nk -i](Xk)) ■ 



The sequence S w [ nk -i](xk) has a convergent subsequence by Lemma 10, and the limit is 
in A w . Thus, x G Sj (Aw). Lemma is proved. 

Lemma 13. Aw C A a ( w ). 

Proof. Again, if x G A w , then x = lim S w [ nk ](xk)- Clearly, 

(x k ) = Ser(w)[n k -1] (S w (p)(x k )) ■ 

The sequence (S w ( )(xk)) is bounded and a{w)[rik — 1] — > a(w). Lemma is proved. 
Corollary 14. If the string w is periodic, then A w = A a i w ). 
The union of individual attractors A w lies inside of K, 

\JA W CK. (6) 

There are many important cases when this union equals K. 
Lemma 15. We have \J A w = K in each of the following cases: 

a) Operators {Sj} are eventually strict contractions, i.e., there exist a < 7 < 1 and 
an integer M > 1 such that for any finite word w* of length > M the operator S w * 
is a contraction with factor 7. (This condition is automatically satisfied if each Sj 
is a strict contraction.) 

b) S~\K)^K forj = 0,...,N-l. 

c) Each operator Sj is invertible on K. 
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Proof. The inclusion (|6]) is obvious. To prove the equality in the special ) and 

b), pick an x G K. There exists a sequence of points {x k } C K, and a sequence w nk 
of lengths n k increasing to infinity such that x = lim S Wn (x k ). We claim that x G A u , 

fc^oo k 

where u = w ni .w n , 2 . . . w nk . . .. Denote u[m k ] = w ni .w n2 . . . w nk . The lengths of the words 
u[mk] go to infinity. 

In the case a), for every k and any y G K we have 

(v)) 

- d(S Wnk (x k ), S Wnk (z k )) 

where z k = S^^y). Then, d(S Wnk (x k ), S Wnk (z k )) < ^ h d(x k) z k ) < 7'* diam(K), where 
l k is the round down of n k /M. Therefore, d(S Wn (x k ), S u [ mk ](y)) — > 0, as k — > 00. Since, 
lim S u [ mk ](y) G A u , and lim ^r^i^) = lim S w (x k ) = x, it follows that x G Au and 

the inclusion if C |J is proved. 

toes 

In the second case, since Sj 1 ^) D K, for every y E K there exist Zj G if with 
y = Sj(zj), j G {0,1,..., N — 1}. Therefore, for every k, we can find y k & K such 
that S , u [ TOfc _ 1 ](y fc ) = x k . Then, (j/ fc ) = S^S^^yk) = S Wnk (x k ). It follows that 

Finally, c) is a special case of 6J. This concludes the proof. 

Remark 16. The case c) may seem too restrictive. However, there are many situations 
where the operators Sj are not invertible on X but are invertible on the attractor K . This 
was first observed by Ladyzhenskaya in the case of Navier- Stokes equations, fTCty . The fact 
is due to the invariance of K and, what is called, backward uniqueness property of certain 
parabolic-like equations. 



Although K equals the union of individual attractors in many cases, there are situ- 
ations when K is strictly larger than that union. This is what we call a Gestalt effect. 



This is a new phenomenon. As we have shown in Lemma 15, the Gestalt effect cannot 
occur when operators Sj are contractions. 

Example of a Gestalt effect. 

In this example the state space X will be the space E 2 of one-sided infinite strings of 
0's and l's. There will be two operators, So and Si, defined as follows: 

S (v)=v(2).v, S 1 (v)=v(l).v 

for all v = v(0)v (l)v(2)v (3) • • • G X. The conditions of Theorem [l] are satisfied, so let 
DJl = K x S2 be the global compact attractor of the corresponding dynamics with choice. 
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[Note that the global compact attractor of the system generated by So is the set of all 
strings with period 3, and the attractor of the system generated by Si is the set of all 
strings with period 2.] 

We claim that the sequence u = 000100 is in K but not in A w for any w G E 2 . 
Let v = 001. cr 3 (t>), i.e., the first three symbols of v are 001, and let Wk = 000. ..0001 
with 3k zeros before 1. Then, for every k, S Wk (v) = 0.001001. ..001 with 001 repeating 
k times. Therefore, S Wk (v) — * u as k —>■ oo, i.e., u G K. To show that u does not 
belong to the union [J A w , we argue by contradiction. If u G A s , then there exists a 

sequence G E 2 such that lim S a [n h ](vk) = u, where n k f oo. Therefore, we can find 

k— >oo 

I, such that S s[ni] (vi), S s[ni+l] (v l+1 ), S s[ni+s] (v l+8 ), all begin with 000100100.... Since 
S s \ ni+1 ](vi + i) = S s r ni+1 y..S s \ ni ](vi + i) = 0001001001..., and the action of operators So and 
Si depends only on the first three symbols in the strings, it follows that vi[3] ^ t^ + i[3], 
because if vi[3] = Uj+i[3], then S s [ ni+1 ](vi + i) starts with at least 4 zeros, i.e., 0000100100..., 
which is impossible. Similarly, Uj+&[3] ^ for j, k = 0, . . . , 8, j k. But there can 

be only 8 different three-letter words in 2 symbols. A contradiction. Hence, u does not 
belong to the [j A w . 



2.4 Dynamics with restricted choice 

As in section [2.3. 1[ E denotes the space of one-sided infinite strings on iV symbols, and d^, 
is the metric on E. Let A be a subshift of E, i.e., A is a closed subset of E and a (A) = A. 
Dynamics with restricted choice is defined on the space Xa = X x A by the operator 
& : (x,w) t— > (S w (o\(x), cr(u>)), where the strings w are now taken from A only. 

We assume that X and So, ... , SV-i satisfy our Assumptions 1 and 2. The existence 
of the global compact attractor, 071a, then follows from the abstract result, Theorem 
[5| The assertions 2 and 3 of Theorem [3] are among the general properties of global 
compact attractors, see Theorem |6j Denote by K\ the projection of 971a onto the X 
component. Clearly, is compact. Also, K\ is a subset of the slice K corresponding 
to the full shift E, as in Theorem [2j Because of the invariance property of 971a, for 
every point y G there is a j, one of the symbols 0, . . . , N — 1, and a point x G K\ 
such that y = Sj(x). Define the sets Aj = {x G K\ : Sj(x) G K\}. It is easy to see 
that each Aj is compact and K\ = Aq U Ai U ■ ■ ■ U Aw-i- By construction, we have 
A U Ai U • ■ ■ U A N - X = S (A ) U S^Ai) U ■ ■ ■ U S N -i(A N ^). 

To analyze the slices Ai\(s) = {i£l : (x, s) G 971a}, we follow the argument of the 



corresponding part of section 2.3.3 



Every point (x, s) G 971a is the limit of the form 



0,s)= lim (S Wk (x nk ),a nk (s nk )) 
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where ( bounded sequence in X, (s n ) is a bounded sequence in A, and Wk is the 

prefix of s nk , s nk = uik-cr nk {s nk ). Because 971a is invariant under 6 and we know that the 
unrestricted dynamics has the global compact attractor 271 = K x £, the sequence (x n ) 
can be taken from the compact K, and we may assume that x nk — > G _K\ Also, we may 
assume that the words Wk converge (to some infinite string G A). The strings a nk (s nk ) 
converge to s. Consider all strings u G A such that Wk-u is a string in A for infinitely 
many k. For every such u we will have x G M\(u). 

We see that the number of different slices of the attractor 971a m &y depend on the 
sequence x nk , but more importantly, it depends on what strings can be attached to con- 
vergent sequences of finite words in A. 

With every sequence (wk) of finite words in A we associate the set s((wk)) of one-sided 
infinite strings u G A such that Wk r u G A for some subsequence Wk v In order to prove 
the third assertion of Theorem [3] we will show that, if A is a sofic shift, the number of 
different sets among all s((wk)) is finite. The argument will be similar to the proof of 
Theorem 3.2.10 in [21]. 

Recall that A is a sofic shift if it has a presentation by a finite labeled graph, see [21] . 
This means that there is a directed graph, G = (V,E), with a finite number of vertices, 
V, and edges, E; the edges are labeled by the symbols 0, 1, . . . , N — 1; from every vertex 
begins at least one infinite directed path; the labels of the edges in the infinite directed 
paths form infinite one-sided strings that exhaust exactly all strings in A. 

Lemma 17. If A is a one-sided sofic subshift of S ; then the number of different sets 
among all s((wk)) is finite. 

Proof. Let G = (V, E) be a labeled graph presenting A. Let (wk) be a sequence of 
finite words allowed in A. For each word Wk pick a finite directed path in G presenting 
it. We can find a subsequence, (wk t ), such that all the words Wk e have the same terminal 
vertex in their presentation. If T is such vertex, then uik r u G A for all infinite paths u 
starting at T. Because the number of vertices is finite, we are done. □ 



Remark 18. Even if the number of different sets among all s((wk)) is > 1, the attractor 
971a may be a product, 97Ta = K\ x A, with the same slice for every string in A. 

Indeed, let = 2 and let A consist of the periodic string u = 100100 . . . and its shifts 
a(u) = 00100 . . . and o~ 2 {u) = 0100 .... If (w^) consists of words ending in 00, then the 
only string that can be attached to Wk is u. If (wk) consists of words ending in 1, then 
the only string is o~(u), and for words ending in 10 the only string is a 2 (u). Thus, we have 
three different sets of the form s((wk))- At the same time, the individual attractors A u , 



A a r u ), and A a 2^, are all equal, as we argue in Corollary 14 
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animated 




Figure 5: The golden+even shift and its animation 

One may ask whether 9JTa is always a product. The answer is no, as the following 
example shows. Let A be the intersection of the one-sided golden mean shift with the 
even shift. In other words, A consists of all sequences of Os and Is such that between any 
two Is there are two or a larger even number of Os. A graph presenting A is given on 
Figure [5j We will animate this graph to define the dynamics. First, identify the nodes 
with three distinct points A, B, and C in M 2 , see Figure [H] left, and define X = {A, B, C}. 
Second, define the maps So and Si acting on points as shown by the directed edges labeled 
correspondingly; for example, S (A) = B, Si(A) = A, and S (C) = B. 

Now consider the set A + of non-empty finite words (blocks) of A. We divide A + in 
three classes and correspondingly divide the strings in A into three classes. The first class 
of words in A + consists of the words ending in 1. Such words can serve as prefixes of 
strings starting with an even (or infinite) number of Os. Denote these classes by A^ and 
aA. The second class of finite words consists of the words ending in odd number of Os. 
The strings for which such words can serve as prefixes are the strings starting with an 
odd number of Os. These classes are denoted by A^ and #A. The last class in A + consists 
of words ending in even number of Os. The corresponding strings are those starting with 
1 or with an even number of Os. These are denoted by Aj and cA. By looking at the 
picture of the animated shift, it is easy to identify the possible limits of sequences S Wk (xk) 
when Wk belong to a particular class, while G {A, B, C}. We see that if Wk G A^, then 
the limit set is {A, B}. If Wk G A#, then the limit set is {B,C}. Finally, if Wk G Aj, 
then the limit set is again {B,C}. Thus, there are two different slices in the attractor 
nJt A . One slice is {A, B}, and the other is {B, C}. We have Ma{u) = {A, B} if u G aA, 
and M-^iu) = {B,C} if u G #A u C A. The global attractor Wl\ is a union of the sets 
{A, B} x A A, {B,C} x B A, and {B, C} x c A. 

Another example of different slices appears in numerical results reported in the next 
section. 
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3 Example 



The simplest mathematical model of malaria transmission goes back to Ross and Mac- 
donald. The state of the human-mosquito interaction system is described by the portion 
of infected humans, x, and the portion of infected mosquitoes, y. The change in time is 
described by the following simple system of ordinary differential equations: 

x = ay (1 — x) — r x , , 

y = bx(l — y) — my 

The nature of the positive coefficients a, b, r, and m is discussed in [27]. In particular, the 
coefficients a and b are proportional to the biting rate and the transmission efficiencies 
(infected human to mosquito and infected mosquito to human), r is the recovery rate (in 
humans), and 1/m is the average mosquito life-span. In practice, it is hard to measure 
these parameters. Also, there are many factors that affect their values, see [27], page 8, 
and the values may change in time. 

The state space for the model is the closed square X = {(x, y) : < x < 1, < 
y < 1}. For initial conditions in X the solution stays in X for all t. If the quantity 
Ro — — is < 1, all trajectories starting in X converge to the origin, and the global 
compact attractor consists of a single point, Pi = (0,0). If Po > 1, the equilibrium 
Pi becomes unstable and there emerges the second fixed point, P 2 = (a;*, y*), inside the 
square X, 

ab — rm ab — rm 

x * = ui \ \ ' y * = ~7ir l — v • ( 8 ) 

b(a + r) a[b + m) 

This second equilibrium is stable, and the global compact attractor of the system consists 
of the two equilibria, Pi and P 2 , and of the heteroclinic trajectory connecting them (and 
staying entirely inside X). The number Po, known as the basic reproductive number, 
detects the emergence of epidemics: when P > 1 there is a stable portion of infected 
population. 

We consider a discrete version of equations ([7]): 

x(t + At) = x(t) + At(ay(t)(l-x(t))-rx(t)) () 
y(t + At) = y{t) + At (bx(t) (1 - y(t)) - my(t)) . U 

The time step map (x(t),y(t)) i— > (x(t + At),y(t + At)) maps X into itself provided 

At < min{^— , — *— } . (10) 

The fixed points for ^ are the same as for ([7]). As in the continuous case, if ab > rm and 
the time step satisfies (10), the global attractor for ^ consists of the two fixed points, 



Pi and P 2 , and the heteroclinic trajectory connecting them. 
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Figure 6: Three individual attrac- 
tors Aw for At = 0.05: left: w = 
111...; middle two: w = 1010...; 
right: w = 000.... 
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Figure 7: Three individual attrac- 
tors Aw for At = 0.005: left: w = 
111...; middle two (very close to- 
gether): w = 1010...; right: w = 
000.... 



We choose two sets of parameters, pset = {a = 4, b = 6, r = 1, m — 2} and 
psetj^ = {a = 2, b = 10, r = 3, m = 2}, and denote the corresponding time step maps by 
5*o and Si. These sets of parameters are not related to any real-life situation but rather 
chosen to better visualize the attractors. The fixed point P2 for pset is (11/15, 11/16) and 
for pset x it is (7/25, 7/12). Figures[2]through|9]show the results of numerical computation. 
The results depend on the size of the time step At. On figures [6] and [7) the left line (the 
heteroclinic trajectory) is the (global compact) attractor for the discrete system (X, Si), 
and the right line is the attractor of (X, So). The two lines between them form the 
individual attractor Aw corresponding to the periodic string w = 1010... (on figure [7] 
the two line are very close). For our example of dynamics with choice, E is the space 
of one-sided infinite strings of symbols and 1. According to Theorem |2j the global 
compact attractor for (X, E) has one slice, i.e., DJl = K x E. The set K for At = .05 
and for At = .005 are depicted on figures [6] and [7], respectively. We have also looked at 
the dynamical systems corresponding to convex combinations of the parameter sets pset 
and pset l and plotted their global attractors. The result is different from K, see figure [6] 
where the "convex combination" is superimposed onto the set K. 

When At — > 0, the upper part of the boundary of K becomes smooth. Note that the 
limit set is not an attractor of any system (|9]) with a fixed, averaged set of parameters 
a, b, r and m. It would be interesting to understand whether the limit set can be obtained 
as a union of the attractors of the systems (X, St), where the operator St corresponds to 
a certain parameter set pset t for some curve connecting pset with pset x in the space of 
parameters. 

Next, we consider restricted dynamics associated with the golden mean subshift A 
(made of one-sided strings of 0s and Is such that each 1 is necessarily followed by 0). The 



graph representing the golden mean shift is shown on figure 10 
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Figure 8: "Convex combination" 
superimposed over K; At = 0.005. 



Figure 9: Golden mean, full; At 
0.05. 




Figure 10: The golden mean shift. 



Our analysis in section 2A shows that the global attractor of the restricted dynamics, 
(X, A) may have at most two different slices: one corresponding to sequences of words 
ending in 1 (the red slice), and the other one corresponding to sequences of words ending 
in (the blue slice). Our computation shows that the attractor of the restricted dynamics 
(X, A) indeed has two slices. The slices are shown on figures 



11 



and 12 



As point sets on the plane, the slices overlap. Their union is plotted on figure M 
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Figure 11: The red slice; At = 0.05. 



Figure 12: The blue slice; At 
0.05. 
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